Spontaneous symmetry restoration in a field theory at finite chemical potential in a toroidal topology 
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We consider the massive vector ^-component (\(p 4 )n theory defined on a Euclidean space with a toroidal 
topology. Using recently developed methods to perform a compactification of a d-dimensional subspace at finite 
chemical potential, we treat jointly the effects of temperature and spatial boundaries, setting forth grounds for an 
analysis of spontaneous symmetry restoration driven by temperature and spatial boundaries as a function of the 
chemical potential. We restrict ourselves to d = 2, which corresponds to the heated system confined between 
two parallel planes (separation L) in dimensions D = 3 and D = 4. We present results, in the large- N limit, 
which exhibit how finite size and chemical potential affect spontaneous symmetry restoration. 

PACS numbers: 11.30.Qc; ll.10.Wx; ll.10.Kk 



I. INTRODUCTION 

We investigate spontaneous symmetry restoration induced 
by both temperature and spatial boundaries. In particular, we 
are interested in studying the influence of a finite chemical 
potential in the context of finite-size effects. We shall fo- 
cus on the A-component vector ip 4 model within the large- A 
approximation for scalar fields. The symmetry restoration is 
carried out by starting from the two-particle irreducible (2PI) 
formalism HI 01 in the Hartree-Fock approximation, which 
sums up all contributions coming from daisy and superdaisy 
diagrams. We need then to implement the 2PI formalism in a 
space with compactified dimensions at large A, which allows 
to obtain nonperturbative corrections to the coupling constant. 
For the sake of simplicity, we study the particular case of fi- 
nite temperature and one compactified spatial coordinate. We 
can then, starting from the broken symmetry region, show the 
behavior of the renormalized mass for different values of the 
temperature, compactified spatial dimension and chemical po- 
tential. 

Field theories defined on spaces with some, or all, of its 
dimensions compactified is of interest in several branches of 
theoretical physics. For a Euclidean D-dimensional space, 
this means that its topology is of the type Tj-, = (S ,1 ) d x 
]g>-D-<^ or j ts coun terpart with the Minkowski signature, with 
1 < d < D, d being the number of compactified dimensions. 
Each of these compactified dimensions has the topology of a 
circle S 1 . We refer to as a toroidal topology. These theo- 
ries are often associated to extra spatial dimensions, as in par- 
ticle physics, where theories defined on toroidal spaces with 
extra spatial dimensions are employed, for instance, as a way 
to investigate the electroweak transition and baryogenesis |3(- 
8]. Also, recent works involving the idea of extra spatial di- 
mensions have been performed in low-energy physics Hd. 

On the other side, in many cases, one is concerned with the- 



ories defined on a space with a compactified subspace. A very 
important development of this kind, which has its roots in the 
late fifties, is the systematic approach to quantum field theory 
at finite temperature, as proposed in 111 lMl4tl . An analogous 
formalism can be constructed for compactified spatial coordi- 
nates, in a D-dimensional Euclidean space. This is an idea 
advanced for instance in 11511 . In this case, as it is stressed 
in 11611 . the compactification of spatial coordinates can de- 
scribe systems confined to limited regions of space, which are 
interpreted as representing samples of material in the forms 
of films, wires or grains. In the general case of quantum field 
theories in toroidal topologies, that is, with compactification 
of the time coordinate and of spatial dimensions, mathemat- 
ical basis to deal with this situation on general grounds are 
consolidated in recent developments Il7lll8ll . This provides a 
general framework for the results from previous works where 
systems at finite temperature and/or compactified spatial di- 
mensions were considered, both for bosonic s yste ms, as for 
instance in Jl9l - l22ll and for fermionic ones j23Tl27ll . 

The paper is organized as follows. In Section 2, symmetry 
restoration is investigated along the lines of 10], which estab- 
lishes a formula for the renormalized mass. This formula is 
then rewritten to take into account the compactification of the 
imaginary time and one of the spatial dimensions, consider- 
ing the dependence with the chemical potential. We can then 
study the behavior of the renormalized mass starting from the 
broken symmetry region. In Section 3, corrections to the cou- 
pling constant, due to compactification as well as the presence 
of a chemical potential are discussed. In Section 4, we col- 
lect the results from the previous sections and, specializing to 
spacetime dimensions D = 3 and 4, we fully determine the 
behavior of the critical temperature as a function of the size of 
the system and the chemical potential. We also present, in the 
same section, our conclusions and final remarks. 
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II. SYMMETRY RESTORATION IN A TOROIDAL SPACE 

We consider the model described by the Lagrangian density 
1 1 u 

C = -d^tfad^fa + ^m 2 (f a (p a + — ((pa^Pa) 2 (1) 

in Euclidean Z?-dimensional spacetime, where m and u are 
respectively the zero-temperature mass and the coupling con- 
stant in absence of boundaries and at zero chemical potential. 
Also, in Eq. ([T), cp a are scalar fields, the sum over repeated in- 
dices is carried from a = 1 to TV and v = 1 to D. We consider 
the large-iV regime, such that linijv->ao mo (Nu) = A, with 
A fixed. To simplify the notation, we drop out the a indices, 
summation over them being understood in field products. In 
order to approach symmetry restoration for this model we first 
follow the 2PI formalism developed in JH @]. In this case, 
one finds a stationary condition for the effective action in the 
Hartree-Fock approximation which translates itself into a gap 
equation, 

G-\x,y) = D-\x,y) + | G(x, x)5 A {x - y), (2) 

where the Fourier-transformed propagators, D(k) and G(k), 
are given by 



D(k) 



1 



k 2 



G(k) 



1 



k 2 + M 2 ' 



(3) 



Here, = (0 \ <p\ 0) is the expectation value of the quantum 
field ip and M is a momentum-independent effective mass. 

In the 2PI formalism, the gap equation corresponds to the 
stationary condition and as such the effective mass depends 
on 4> and conveys all daisy and superdaisy graphs contributing 

toG(fc) una. Nevertheless, in order to investigate symmetry 

restoration, we can take instead a particular constant value M 
in the spontaneously broken phase. Renormalization of the 
mass m and of the coupling constant u can be performed along 
lines similar to those in 1I2T I28I1 . leading to the equation 



M 2 = -m\ 



2 V 



G(M), 



(4) 



where m 2 R and ur are respectively the squared renormal- 
ized mass and the renormalized coupling constant, both at 
zero temperature and zero chemical potential, in the absence 
of boundaries, and G(M) is the finite part of the integral 
G(x,x) = (1/(2tt) d ) / d D kG(k), which will be calcu- 
lated using dimensional and zeta-function regularization tech- 



niques. Notice also the minus sign of the m 2 R term, which is 



the same choice made in (0, H^] to ensure spontaneous sym- 
metry breaking. Therefore, Eq. (|4]i will be our starting point. 
It gives the value of the renormalized mass at the broken- 
symmetry phase and can be rewritten as 



m 2 ((j)) 



-M + ^ 



Ur 

2 



d D k 



(2ir) D k 2 +M 2 ' 



(5) 



where the effective renormalized mass fh 2 ((j)) — —m R + 
(ur/2)4> 2 has been introduced. In the sequel we will obtain 
the generalization of the above equation in such a way as to 
include the toroidal topology as well as the chemical potential, 
but first we notice that restoration of the symmetry will occur 
at the set of points in the toroidal space where m 2 is null. 

We now proceed to generalize Eq. ® to a theory defined on 
a space with a toroidal topology. In the general case, the sys- 
tem is in thermal equilibrium with a reservoir at temperature 
and confined to a [d — 1) -dimensional spatial rectangular 
box of sides Lj, j = 2, 3, d. We use Cartesian coordinates 
r = (xi, Xd 1 z), where z is a (D — d) -dimensional vector, 
with corresponding momentum k = (fci, fed,q), q being 
a (D — d) -dimensional vector in momentum space. Then the 
Feynman rules should be modified according to fHHHI 



dk T 
~27 



1 OO 

5 E 





2^z- £ 



n T — — 00 
+ 00 



2n T 7r 
k T — > — ifj,, 



2n.;7T 

Kj > — - — , l = 2, o, u, 



rii — — 00 



(6) 



where r corresponds to imaginary time and fi is the chemical 
potential. We consider the simpler situation of the system at 
temperature /3 _1 and one compactified spatial coordinate (^2) 
with a compactification length L2 = L. In this case, using 
Eq. ©, we can perform a suitable generalization of the proce- 
dure in J2l, to take into account finite-size, thermal and bound- 
ary effects in Eq. (f5]). The integral over the £>-dimensional 
momentum in Eq. (O becomes a double sum over n T and 
n-2 = n x together with a (D — 2)-dimensional momentum 
integral. Then, following steps similar as in 1U9I1 and using 
dimensional regularization to perform the integral l29tl . the 
renormalized ({3, L, /i)-dependent mass in the large-iV limit 
can be written in the form 



m 2 (T,L,ii) = -M 



2 , Ur 



M D-2 n (D-2)/2 r ( s _ D_l) 



4tt 2 r (s) 



-\J a T a x 



E 



T I ' b T 





—li ) + a x n x + c 



(D-2)/2-s 



(7) 
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where we have changed variables in the integral, ki/2ixM — >• 
q,i, and introduced the dimensionless quantities a T = 
(M/3)" 2 , a x = (ML)- 2 and c = (2tt) The double 
sum in Eq. (0 is recognized as one of the inhomogeneous 



Epstein-Hurwitz zeta functions, Z 2 (s ~ -^j -2- ; a r .a x ; b T , b x ), 
which has an analytical extension to the whole complex s- 
plane J30|,[3ll]; m general for j = 1, 2, 



Z c 2 (v;{a J };{b j }) = 



tt\c\ 2 - 2u T(v- 1) 47T V |C| 1 



j=l n 3 =l 



j 2lTCTlj 

V V«7 



2 2J cos(27mi6i)cos(27m 2 62) \f+~ 

n i ,n2 — — oo 



V ai a 2 



(8) 



For us, a\ = a T , a 2 = a x , b\ = b T = i/3fi/27t, b 2 = b x = 0, 
c = 1/2tt and u = s — (D — 2)/2. Replacing Eq. © into 
Eq. (0, the thermal and boundary corrected mass is obtained 
in terms of the original variables, f3, L, li and of the fixed 
renormalized zero-temperature coupling constant in absence 
of boundaries, 

X R = lim (Nu R ). 

However, the first term in Eq. ([8]l implies that the first term 
in the corrected mass is proportional to r(l — D/2), which is 
divergent for even dimensions D > 2 IU9I1 . This term is sup- 
pressed by a minimal subtraction, leading to a finite effective 
renormalized mass; for the sake of uniformity, this polar term 
is also subtracted for odd dimensions, where no singularity 
exists, corresponding to a finite renormalization. 

In dimension D, the renormalized zero-temperature cou- 
pling constant in absence of boundaries Xr has dimension of 



mass 4 ~ D ; accordingly, we introduce a dimensionless cou- 
pling constant A'^ = Xr M d ~ a . Also, we define the dimen- 
sionless reduced temperature t, reduced chemical potential tu, 
and the reduced inverse length of the system x, in such a way 
that we have for any dimension D, the set of dimensionless 
parameters defined by 

X' H = X R AI d - 4 , t = T/M , x = L^/M , u = fx/M. 

(9) 

We then can obtain, after subtraction of the polar term, 
which does not depend on /3, L and fx, the corrected mass, 
fh 2 (D, (5, L, fx). This implies that the condition for symme- 
try restoration, fh 2 (D, (5, L, li) — 0, can be written in terms 
of the above dimensionless parameters, replacing A^, by the 
corrected coupling constant X' H (D, f3, L, li) (this is precisely 
defined in the next section, Eq. ( fT5] l) , in such a way that the 
critical equation reads 



{2n) D / 2 




n,l=l 



^(t)+E(?)' 



K D 



Kd 




(10) 



III. CORRECTIONS TO THE COUPLING CONSTANT 



four external legs, which leads to the expression (we consis- 
tently define u' R = u R M D ~ 4 ) iH 



In this section, we follow the reasoning made in 112 1H . ap- 
propriately modified to incorporate the effects from the chem- 
ical potential. We consider the zero-external-momenta four- 
point function, which is the basic object for our definition of 
the renormalized coupling constant. At leading order in 1 /N 
it is given by the sum of all chains of one-loop diagrams with 



i + Nu' H n.(D,p,L,fiy 



(ii) 
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where the dimensionless one-loop diagram is given by 

I— — — oo 

where 

-E-n T n x (^) " 



s=2 



Then, proceeding along lines similar to those leading to 
Eq. ( [Tol l, we write IT(Z?, (3, L, fj,) in the form 



U(D,p,L,n) = H(D) 



(2n) D / 2 



R(D,p,L,n), 



(13) where R(D, /3, L, p,) is given by 




and H(D) oc T (2 — -j) is a polar parcel coming from the 
first term in the analytic extension of the zeta function in 
Eq. ©. Notice that from general properties of Bessel func- 
tions, the above equations are meaningful for a reduced chem- 
ical potential satisfying the condition < oj < 1. This 
is the allowed range of u in Eqs. ( TTOb and (TBI) . We see 
that for even dimensions D > 4, H(D) is divergent, due 
to the pole of the gamma function. Accordingly, this term 
must be subtracted to give the renormalized single bubble 
function Hr[D, j3, L, p). We get, simply, U R (D, p, L, p) = 
[l/{2it) D/2 )]R(D, (3, L, p). As mentioned before, the term 
H(D) is also subtracted in the case of odd dimensions 
D, where no poles are present, corresponding to a finite 
renormalization. Using properties of the Bessel functions, 
we see that, for any dimension D and finite values of the 
chemical potential p, R(D, j3, L, p) satisfies the conditions 
lim/^L-^oo R(D, /3, L, p) = 0, lim^^^o R(D,j3,L,p) -> 
00, and R(D, f3, L,p) > for finite \i and for any val- 
ues of D, (3 and L. Under these conditions, we define 
the dimensionless j3-, L- and /i-dependent renormalized cou- 
pling constant X' R (D, f3, L, /1) at the leading order in 1/N as 

X' R {D,/3,L,fi) = NT { n\Q,P,L,ii), which, from Eq. O, 
after subtraction of the polar term H(D) in H(D, j3, L, 
leads to 

X> R (D, f3,L,(j,) = — Vfi[1/(2)r)fl/ a Jp(A ^ ^ ~ } ■ d5) 

In the next section we will investigate the restoration of sym- 
metry, taking into account thermal, boundary and finite chem- 
ical potential corrections to the coupling constant as presented 
above. 



IV. BOUNDARY AND CHEMICAL POTENTIAL EFFECTS 
ON THE SYMMETRY RESTORATION: COMMENTS AND 
CONCLUDING REMARKS 

In the general situation, Eq. ( TTOb does not allow an alge- 
braic solution. For numerical evaluations, we fix the value 
A^j = 0.50 and take several values of the dimensionless pa- 
rameters t, \ and uj. We present and comment on the results 
conccurently for dimensions D = 3 and D = 4, since the 
behavior of the system is qualitatively the same, apart from 
numerical details, in both cases. 

In Figs. Q] and [2] we exhibit the critical temperature as a 
function of the reduced inverse size of the system for different 
values of the chemical potential, for D — 3 and D = 4 re- 
spectively. We see from them that the behavior of the critical 
temperature is different, by changing the values of the chemi- 
cal potential, for small and large values of \ (large and small 
sizes of the system, respectively). An interesting aspect, ex- 
plicitly shown in Figs. Q] and [2] is the existence of a particular 
size of the system, Lo, corresponding to the reduced inverse 
size xo> where the critical temperature vanishes. This particu- 
lar value xo is independent of the chemical potential for both 
D = 3 and D = 4. This is emphasized in the lower plots 
of Figs. Q] and [2] which show in detail the domain around the 
characteristic value x = Xo- F° r me value of the reduced cou- 
pling constant we take (A^. = 0.50), we find xo ~ 69.96 for 
D = 3 and xo ~ 7.78 for D = 4. 

Although these results seems a priori unexpected, we can 
find the above values of Xo> f° r b° m D = 3 and D = 4, 
directly from Eq. ( fTOb . Indeed, it should be noted that for 
X = Xo< tri e symmetry breaking region disapears completely 
and we have a null critical temperature. Then xo is obtained 
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FIG. 1: Reduced critical temperature as a function of the reduced 
inverse size of the system for dimension D = 3 (upper plot). We 
fix X' R — 0.5 and take the chemical potential values ui = 0.1 (full 
line), 0.3 (dashed line) and 0.4 (dotted line). The symmetry-breaking 
regions are in the "inner" side of each curve. The lower plot is a 
"zoom" enhancing the region of the characteristic size of the system, 
corresponding to \ ~ 69.96. 

by solving Eq. ( fTOb for t = 0, in which case only survives the 
second term, in such a way that all dependency coming from 
the chemical potential drops out. Therefore, for this particular 
point, t = 0, for all sizes lower than a characteristic value, 
the chemical potential magnitude does not have any influence. 
Essentially what appears to happen, is that in this case, the 
finite-size behavior of the physical system collapses into the 
one corresponding to a zero chemical potential, as is the case 
for a Bose-Einstein transition. 

We also see from Figs Q] and [2] that, for each value of ui, 
there is a limiting smallest size of the system, L m j n (w), corre- 
sponding to a largest reduced inverse size x max (w), such that 
Xmax(w) > Xo> over which the transition ceases to exist. 

Moreover, we can see clearly from Figs.[T]and[2]that, by ef- 
fect of the spatial boundaries, \o is the border between two re- 
gions: x < Xo an d Xo < X < Xmax, with different behaviors. 
In the first region, the critical temperature is uniquely defined 
in terms of the size of the system and of the chemical poten- 
tial: For each pair (x, ui) there is only one critical temperature, 
while in the second region, two values of t c may exist for the 
same values of u> and x- I n tne region Xo < X < Xmax, 



FIG. 2: Reduced critical temperature as a function of the reduced 
inverse size of the system for dimension D = 4 (upper plot). We 
fix the (dimensionless) coupling constant Xr = 0.5 and take the 
chemical potential values ui — 0.0 (full line), ui = 0.5 (dashed line) 
and oj — 0.9 (dotted line). The symmetry-breaking regions are on 
the "inner" side of each curve. The lower plot shows in detail the 
critical temperature in the region around the characteristic size of the 
system. 

there are for each value of ui, two possible critical temper- 
atures, say, tc an d t c , with tc > t^\ associated re- 
spectively to the lower and the upper branches of the critical 
curve. This means that in this region we have two possible 
transitions. For D = 3, we take, from inspection of Fig.Q] 
t = ti^ = 1.93, t = = 6.89 and an intermediate temper- 
ature, tc < t < t c 2 \ t = 4.80. For these temperatures we 
plot in Fig. [3l curves of the effective potential of the system, 
given by U(J3, L, fi) = ±m 2 (A L, ^ 2 + ±X' R (j3, L, ^ 
(we define the quantity <p by taking ip 2 — (fia'fia), for fixed 
values of the reduced chemical potential, us = 0.30 and of 
the reduced inverse size of the system, x = 100. These plots 
confirm that two of the temperatures, t = ti = 1-93 (full 

(2) 

line) and t = t c = 6.89 (dashed line), are critical tem- 
peratures corresponding respectively to the lower and upper 
branches of the critical curve; the intermediate temperature, 
t c < t < t c 2 \ t = 4.80 (dotted line), corresponds to the sys- 
tem in the symmetry-broken region. A similar kind of "dou- 
bling" occurs for large sizes (but not for L —> oo) of the sys- 
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FIG. 3: Curves for the effective potential U = -|m 2 </? 2 + -ii^'af 4 , 
for fixed values of the reduced chemical potential, u> — 0.30 and 
of the reduced inverse size of the system, \ — 100 for D — 3; ip 
and U are measured in units of M 5 and M 3 , respectively. We take 
three values of the reduced temperature, two of them being critical 
temperatures, to = 1-93 (full line) and t c 2 ^ = 6.89 (dashed line), 
and an intermediate one, t — 4.80 (dotted line), in the symmetry- 
broken region. 



tern. In this case two different sizes of the system with the 
same chemical potential may correspond to the same critical 
temperature. In the region \ < Xo> the critical temperature 
grows for increasing chemical potential for fixed \ an d a ls° 
for the system in unlimited space (L — > oo). In this case, we 
may find the bulk reduced critical temperature ij? ulk (^) taking 
the limit \ -» in eq. ( TTOt . 

As an overall conclusion, the results suggest that finite-size 
effects with finite chemical potential are relevant and deeply 
changes the critical curves with respect to the ones for the 
system in bulk form. In particular, these actors lead to the ap- 
pearance of a "doubling" of critical parameters, which, up to 
our knowledge, is not a trivially expected behavior. This be- 
havior is to be contrasted with what happens with the system 
in bulk form, where there is always an unique critical tem- 
perature, which grows with increasing chemical potential. We 
also show the existence of a characteristic size of the system, 
which determines the existence, for each value of the chem- 
ical potential, of a minimal size sustaining the broken phase. 
This characteristic size, obtained from Eq. J10I . is the same 
for all values of the chemical potential. 
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